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1. Introduction 

Quantum dimer models (QDMs) are among the simplest lattice models known to display 
exotic quantum phases [U [21 El SI E]. In addition to conventional, symmetry broken 
dimer crystal phases, QDMs have quantum liquid phases with no local order. In 
particular, the QDM on the triangular lattice has been shown to possess a quantum 
liquid ground state with topological order [3l [6] . Unlike conventionally ordered phases, 
topological phases have no conventional broken symmetry and therefore can not be 
distinguished by a local order parameter [71 E]- Features of topological order include 
a robust topological ground state degeneracy and a gap to anyonic quasiparticle 
excitations. Because of these unique features, topological phases have been proposed as 
the basis for a physically fault tolerant quantum computer [9l [TOl [8] . Several spin and 
Bose- Hubbard models are known to map to quantum dimer models [HI [12]; consequently 
there have been proposals to realize triangular lattice QDMs in various experimental 
systems as the basis for a topologically protected qubit [I3l [TH [15] . 

A transparent physical picture of topological order exists for lattice models with 
loop degrees of freedom. Such quantum loop model models are known to possess 
topologically ordered phases [101 [lEl [TTl [TBI [IE] , and the formation of topological order 
can be described in terms of loop condensation. In the loop condensation picture, the 
creation of topological order coincides with the formation of fluctuating loops on all 
length scales. While the microscopic degrees of freedom in QDMs are non-overlapping 
dimers, QDMs can be mapped to closed loop models [20l [21]. In this paper we use the 
mapping of QDMs to closed loop models as a tool to analyze the phases of the triangular 
lattice QDM. In particular we consider the geometric properties of the loop distribution 
to understand the loop nature of the liquid and crystalline phases of the QDM. 

The structure of this paper is the following: first we review the triangular lattice 
QDM and the concept of loop condensation and topological order in sections [2] and [3l 
We then introduce the loop mapping of the quantum dimer model in section [H The 
main results of this paper are then presented: we study several wavefunctions that are 
representative of phases of the triangular lattice QDM and capture transitions from the 
dimer liquid to dimer crystal phases. Using Monte Carlo sampling of these wavefunctions 
with local and directed-loop updates, we investigate the loop nature of these phases 
in sections [5B This allows us to characterize the loop properties of the liquid and 
crystalline phases of the triangular lattice QDM. 

2. The quantum dimer model on the triangular lattice 

First we briefly review the quantum dimer model on the triangular lattice. In a hard- 
core dimer model, the degrees of freedom are dimers that live on the links of a lattice, 
and the hard-core constraint forbids more than a single dimer from touching each vertex. 
The quantum dimer model was first introduced on the square lattice by Rokhsar and 
Kivelson [1], and later generalized to the triangular lattice by Moessner and Sondhi [3]. 
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The Hilbert space of the QDM comprises fully-packed, non-colliding dimer coverings of a 
lattice, where all dimerizations are by definition orthogonal. The only way to rearrange 
dimers in a fully-packed dimerization without violating the hard-core constraint is to flip 
dimers around a closed loop on the links of the lattice, where the links in the loop are 
alternately occupied and unoccupied by dimers. Consequently, the minimal quantum 
dynamics give a resonance between the two orientations of dimers around a "fiippable" 
plaquette. On the triangular lattice a fiippable plaquette is defined as a length 4 
rhombus lattice containing two parallel dimers. The canonical QDM Hamiltonian on 
the triangular lattice [H [3] is: 

i/RK = -^(i^x-^i + ^-c-) + ^(i^x^i + i-^x-^i). (1) 

p 

In ([1]), the sum is over all rhombus plaquettes labeled by p, including all three 
orientations, and \^) and represent the two possible fiippable dimer configurations 
around p. The t term gives the dimers kinetic energy, and the v term represents an 
interaction between parallel dimers. In this work we will only consider t > 0. 

These local dynamics are not ergodic over all dimer configurations, and 
consequently split the Hilbert space into topological sectors that are defined by quasi- 
ergodicity of the local dynamics. This means that on the triangular lattice there are 
a finite number of topological sectors, within which the plaquette flip dynamics are 
ergodic [3]. On a lattice with periodic boundary conditions, one can also classify 
topological sectors by the parity of the number of dimers crossing a closed loop that 
passes through the faces of the plaquettes and encircles the lattice; a plaquette flip 
conserves this parity. On a torus, there are 4 such parity sectors, two for each of the 
directions around the lattice. However, there also exist a finite number of symmetry 
related "staggered" configurations with no fiippable plaquettes, which are frozen under 
plaquette flip dynamics; it is believed that plaquette flips are ergodic in each parity sector 
excluding these staggered configurations [3]. We label these parity winding sectors as 
{E,E), {E,0), {0,E) and (0,0), where E and O represent the even and odd parity 
sector, respectively. 

At the point v = t, the so-called RK point, the exact zero energy ground state of 
-f^RK is the equal superposition of all dimerizations: 

\RK) = J2 \C) (2) 

c 

where the sum in equation (|2]) is taken over all dimer configurations O in a topological 
sector [1]. Since the norm of \RK) is equal to the partition function of the corresponding 
classical dimer model, expectation values of diagonal observables are identical to those 
of the classical dimer model. Dimer-dimer correlation functions are known to decay 
exponentially in the classical dimer model on the triangular lattice [221123] . and so \RK) 
describes a dimer liquid phase. Additionally, imaginary-time correlation functions of 
equation (jj]) can be related to dynamic correlation functions of a Monte Carlo simulation 
of the classical dimer model [211 [251 126]; on the triangular lattice, such calculations 
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have shown that -/^rk has a finite gap at the RK point Quantum Monte Carlo 

calculations by Moessner and Sondhi [3] and Ralko et. al. [6] demonstrate that the 
dimer liquid phase extends beyond the RK point for a finite range oi v/t for v <t and 
possesses Z2 topological order. While the QDM has been related to Z2 gauge theory in 
the literature, no exact mapping from the deconfined phase of a Z2 gauge theory to the 
dimer liquid is known [281 129] . 

For V > t, the triangular lattice QDM is in a staggered dimer crystal phase, as the 
staggered configurations are frozen, zero energy ground states of -ffRK- For l^;! >> t and 
v/t < 0, states with the maximum number of fiippable plaquettes are favored. The 12 
symmetry related columnar configurations (see figure [1]) have the maximum number of 
fiippable plaquettes, Ni/Q, where A^; is the number of links in the lattice. Additionally, 
any configuration reached by translating dimers along any number of rows of a columnar 
configuration and by rotating all dimers in any number of columns of a columnar 
configuration, generates another maximally fiippable dimerization[3]. References [3j and 
[6J show that quantum fiuctuations favor the columnar order, and that for v/t < —0.75, 
i^RK has a columnar dimer crystalline ground state. The staggered and columnar phases 
break both translational and rotational symmetries of the lattice. 

Between the columnar crystal and the dimer liquid lies a resonating dimer crystal 
phase with a 12-site unit cell (see figured]) termed the a/T2 x y/T2 phase [31 [301 E]- In 
this phase, dimers resonate within 24 link hexagons, and most translational symmetries 
are broken. However this phase retains a rotational symmetry of the lattice. Ralko 
et. al. showed that the dimer correlation functions are qualitatively reproduced 
by a wavefunction that is an equal superposition of all dimerizations within the 
hexagons [30l[29|. 




(a) (b) 

Figure 1. Reference dimer configurations that are representative of dimer 
crystals in the triangular lattice QDM. The thick red dimers represent the reference 
configurations, Rq. (a) A columnar configuration that defines Rcoi- (b) In the 
y/12 X VT2 phase dimers resonate in 12 site hexagons, colored red in this figure. The 
thick red dimers show the reference configuration R^/i2- 
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3. Loop condensation 

Here we briefly review the paradigm of loop condensation in topological phases. Loop 
condensation has been discussed as a physical mechanism for the formation of topological 
order by several authors [Ml HOl UHl HTJ. Quantum systems with loop-like degrees of 
freedom may display a variety of phases. These include symmetry broken loop crystal 
phases for which the length of loops is peaked at particular length scale, as well as 
scale-invariant loop liquid phases characterized by fluctuating loops on all length scales. 

One of the simplest exactly solvable models with a quantum liquid ground state is 
the toric code fl0\ : 

H^c = -KY,A,-\^Y,B^, ^ = J]a|, Bp = J]aJ. (3) 

Here, {uj} are spin-1/2 degrees of freedom that are located on the links of a square 
lattice on a torus, v and p label the vertices and plaquettes of the lattice. When 
= the ground state, I^E'tc), is an eigenstate of all A^, i.e. \^tc) = + I^tc) ior 
all V. These eigenstates of the vertex operators A^ have 0, 2 or 4 down spins touching 
each vertex. The ground state subspace can be interpreted in terms of a closed loop 
model by choosing a reference configuration \R) in the basis (from the subspace of 
+1 eigenstates of all A^) that is defined as the empty loop state. For example, we may 
choose \R) to be the spin polarized state, |{V(t| = +1}). Any other configuration in the 
Ay = +1 subspace can be reached by applying a product of closed loop operators: 

\A) = n ^^1^)' ^^^n^i' 

where £ is a closed loop along the links and {^aIr a closed loop covering of the square 
lattice. Since [A^j, We] = 0,Wv, \A) is in the A^ = +1 subspace for all {^a}- In this way, 
the closed loop covering defining \A) is 

In general {iA}ji may comprise both contractible loops and non-contractible loops 
that wind around the torus. The winding sector {wi,W2)a is defined by the number 
of windings in {iA}^ about the two axes of the torus. Though the the plaquette term 
Bp causes fluctuations in these loop coverings, the parity of the winding number is 
conserved. The operator 

Wi,2 = n ^1 (5) 

ieci,2 

measures the parity of the winding sector and commutes with H^c, where ci,2 are loops 
that pass through the faces of the plaquettes and wind about one of the axes of the 
torus. Therefore, the ground state subspace is divided into four topological sectors that 
are defined by the parity of these two winding numbers. 

The ground state |^tc) is an equal superposition of all loop coverings in a given 
topological sector, and has fluctuating loops on all length scales. This quantum loop gas 
has no local order, but does possess topological order and is described by the deconfined 
phase of Z2 gauge theory [lO]. By considering the limit Ag — )■ cxd, perturbations can drive 
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transitions within the closed loop subspace [SI]. For example addition of a magnetic 
field H' = —h'Y^. 0"! will drive a transition to a spin polarized phase: for h > \hc\ the 
ground state is a dilute loop crystal, and for /i < — \ hc\ the ground state is a fully packed 
loop crystal. These transitions do not involve spontaneous symmetry breaking; we can 
alternately consider the effect of adding a loop interaction such as H' = JYl,{ij) ^t^j 
which favors a rotational symmetry broken loop crystal for J ^ A^- We note that 
loop crystals need neither be dilute nor involve short loops; however, we may choose 
\R) to reflect the broken symmetries of the crystal phase, such that the resulting loop 
configuration will involve only short loops. 

4. Loop order in quantum dimer models 

A dimerization of a lattice, |C), can be interpreted as a closed loop covering of the lattice 
by superimposing it with an arbitrary reference dimerization, \Rq) (see flgureE]) [20|l2T]. 
The resulting doubled dimerization, {C : -Rq}; known as the transition graph, has both 
a reference and physical dimer touching each vertex; this therefore forms closed loops 
on the lattice with the exception of where dimers in \C) coincide with those of |-Ro)- It 
is necessary to choose a convention for the deflnition of such coinciding dimers. Two 
possible interpretations of overlapping dimers are: (a) they are considered to be a loops 
of length 2, or (b) they are deflned as an empty links of the transition graph. Choosing 
option (a) means that the transition graph is a fully packed loop covering of the lattice 
and the transition graph {Rq : Rq\ comprises all length 2 loops. Choosing (b) means 
that the transition graph is not fully packed, and that {Rq : Rq] is an empty loop 
covering. In order to make an analogy with the deflnitions of loops given in equation 
dl]), we will generally choose below the convention (b), except where noted. 




Figure 2. A transtion graph loop configuration {C : Rq} formed by superimposing a 
pliysical dimer configuration \C) (black) with the reference dimerization \Rcoi) (red). 

Choosing a reference dimerization |i?o) translates a dimer conflguration into a loop 
conflguration. However, unlike the toric code, there is no reference conflguration that 
retains all of the lattices symmetries. Therefore we are forced to choose an \Rq) that 
breaks translational and rotational symmetries of the lattice. Certain properties of the 
resulting loop model will depend on this choice of a symmetry broken \Rq). All other 
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dimerizations can be reached from \Rq) by applying a product of closed loop operators 
that alternately pass through the links that are occupied and unoccupied by dimers in 
|i?o): 

|C) = H W, \Ro) , W,= II d-dt. (6) 

In equation ([6]), d'j'{dj) are dimer creation (annihilation) operators and (r, s) are a pair 
of links in i that meet at a vertex, that are occupied and unoccupied, respectively, 
by a dimer in |-Ro)- The set {ic}Ro is simply the transition graph {C : Rq}, where 
we interpret hnks that are occupied in both in \C) and \Ro) as empty links in {ino}- 
We note that the choice of |-Ro) restricts the possible loop coverings such that if a loop 
touches a vertex, it must pass through the link occupied in |-Ro)- This dimer-loop Hilbert 
space is not equivalent to the Hilbert space of all loop coverings of the lattice. 

The phases of -f^RK may now be re- interpreted as those of a loop model. \RK) is 
the equal superposition of all loop coverings that pass through the dimers of |-Ro), and 
therefore is a loop gas. We expect the topological phase adjacent to the RK point to 
be described by a scale-invariant loop liquid. Static dimer crystal phases, such as the 
columnar and staggered phases, correspond to loop crystals with a narrow distribution 
of loop lengths; if \Rq) is chosen to be a dimer configuration with the corresponding 
broken symmetry, then the loop model is in a short looped dilute phase. The loop 
order in the resonating a/12 x vT2 phase is less obvious: if \Ro) is chosen to be one of 
the dimerizations with vT2 x \/T2 order, the ideal a/T2 x a/T2 state will be a dense, 
fluctuating short looped phase. This transition loop description of the dimer model 
allows us to make a connection to the loop condensation picture of topological order. 
A dimer crystal to dimer liquid transition is then described as a transition from a loop 
crystal to a scale-invariant loop liquid. 

5. The dimer liquid at the RK point 

Here we analyze the geometric properties of the loop description of the RK wavefunction, 
taking this to be representative of the dimer liquid. This analysis follows a similar 
approach as that of Sutherland and Kohmoto & Shapir in their studies of the spin- 1/2 
RVB wavefunction on the square lattice [201 EH ESI EHl IMl EH] . To understand the loop 
properties of \RK), we can refer to the classical 0(n) loop model [361 EH EH]. In this 
model, which can be solved exactly on the honeycomb lattice, configurations are closed 
loop coverings of the lattice, and the partition function is given by 

c 

where C is a closed loop configuration comprising TV loops of total length C In equation 
(I7j) , is the weight per length of loops and n is the loop fugacity. For n <2, there are 
two phases separated by a critical line at Kc = [2 + v^2 — n]~^^'^. separates a dense 
loop phase for K > Kc and a dilute phase for K < K^.. The dense phase is characterized 
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by a power law distribution of loop lengths, (s) ~ where Pi (s) is the density 
of loops of length s [39]. Additionally, there is a large spanning loop Em with fractal 
dimension Df, whose length scales with the system size L as s(£m) ~ [40]. These 
exponents are related by the scaling dimension x: 

Dj = 2-x, p=l + ^ (8) 

Within the dense phase, x is given by 

a^ = |-^(^-l)% n = -2cos(7r^) (9) 

for < g < 1 where g is the coupling constant of the Coulomb gas description of of the 
0{n) model [38]. The 0{n) loop model on the triangular lattice has been studied by 
Knops et. ai, and is known to possess a critical dense phase of the same universality 
class as that of the honeycomb lattice [H]. 

In the 0(1) model on the honeycomb lattice, Kc = l/"\/3, so = 1 is in the 
dense phase. At i^' = 1, ^o(i) is an equal weighted sum over all loop configurations 
with the critical exponents Df = 7/4 and p = 15/7 |3H1 HO]. The RK wavefunction 
can similarly be described as equal superposition of all transition loop configurations 
within a topological sector, for a given choice of \Ro), and therefore we can relate the 
loop configuration at the RK point to the critical phase of 0(1) model. 

To calculate the loop properties of \RK) we choose a reference configuration 
\Ro), and then compute the distribution of loops in the transition graph, P^{C,s) = 
hf°{C, s)/L'^ where hf°{C, s) is the number of loops of length s in {C : Rq} and L is 
the linear dimension of the lattice. We have computed the expectation value of P^° by 
Monte Carlo sampling of \RK) using local plaquette flip updates on triangular lattices 
of up to L = 192, with A^^ = 3 * links. 




Figure 3. Loop distribution Pi vs. loop length s at the RK point on a i = 120 
triangular lattice using reference configurations \Rcoi) and in the (E, E) winding 

sector. The line shows the fit to a power law decay. 
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(a) (b) 

Figure 4. (a) Longest loop distribution, Pi^j vs. loop length s{£m) at the RK point 
using reference dimerizations \Rcoi) and \R^^), in the (i?, E) winding sector on lattice 
size L = 72. (b) Finite size scaling of the expectation value of the length of the longest 
loop s{£m) at the RK point for \Rcoi) and l-Ryx^): in the {E,E) winding sector of 
lattice sizes up to L = 192. Errorbars are smaller than the symbols shown. The line 
shows the fit to s{£m) oc L^f . 



Here we will consider two choices of reference configurations: \Rcoi), the columnar 
reference state shown by the red links in figure [TK, and \R^), which is the configuration 
illustrated by the thick red links in figure [TJo that is one of the equally weighted 
configurations of the ideal vT2 x a/T2 crystal. In figure [3l the loop distribution Pi{s) is 
plotted for these choices of \Ro)- We see a clear power law over a range of length scales. 
The best fit power law p = 2.12 ± 0.01; this is consistent with the theoretical value for 
the 0(1) model, p = 15/7 ^ 2.14 [381 [39], considering the imperfect power law behavoir 
of Pi on lattices of these sizes. Additionally, we see that is independent of the 

choice of |-Ro) for lengths longer than about 16. We have computed the distribution of 
the longest loop Em, PiAfi^), and this is plotted in figure [1^. To compute Dj at the RK 
point, we analyze the finite size scaling of Im, shown in figure [H (b). The best fit gives 
Df = 1.751 ± .001, which agrees with the known fractal dimension of the 0(1) model, 
Df = 7/4 [38l[40]. 

We can also characterize the loop liquid in terms of the total loop density, 
p = C/Cm, where Cm = is the maximum possible total loop length. The value p = I 
corresponds to the maximum loop density, where no dimers reside on links occupied 
in l-Ro), while p strictly vanishes only for |-Ro)- Consequently, in crystalline phases, p 
depends on the choice of |-Ro)- At the RK point, all links are occupied with probability 
1/6, and the total loop length is twice the average number of dimers on links unoccupied 
in I-Rq)) which is 5/6 x Ni/6. Therefore at the RK point, p = 5/6; this agrees with our 
numerically computed value of p = 0.83333(3) on a L = 120 lattice, computed for \Rcoi) 
and l-R^). The RK wavefunction can be described as a dense loop gas. We note that 
in this loop gas phase, both the total density and geometric properties of loops are 
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independent of the choice of \Ro). 

6. Loop order in dimer crystal phases 

Here will consider the loop order in transitions from the dimer liquid at the RK point 
to two dimer crystals that are present in phase diagram of the cannonical QDM as 
defined in equation ([1]). To analyze these transitions, we will consider a wavefunction 
that interpolates between \RK) and a given dimer crystalline order: 

Here we will consider two possible colorings of the triangular lattice where the red links 
define links that are favored by the crystalline order, to which we give weight 1, and 
Nb is the number of occupied black links, which are given a weight z. In equation ( ITO|) . 
Z = '^(j z'^^^^'"\ which is the partition function of a non- interacting classical dimer 
model where black links are given a fugacity z^. Additionally, (^)) is the zero energy 
ground state of a local, stochastic matrix form Hamiltonian that is generalization 

of HjiK- 

H, = J2 -^(l^X^I + ^-c-) + ^(^^^^ \^){^\ + z-^^^ I^X^l) (11) 
p 

where ANb = Nb{^) - Nb{^). For 2 = 1, (^ = 1)) is the RK wavefunction, and 
and for z = 0, |\& (2; = 0)) is an ideal dimer crystal as defined by the red links of the 
lattice. 

As described in section Ej Hyik displays two crystalline phases for t> < t, in addition 
to the dimer liquid phase. These crystal orders are illustrated in figure [TJ the red links 
of panel (a) show one of 12 symmetry related columnar states, and the red links of 
panel (b) are favored by the a/T2 x a/12 resonating crystalline order. We will use these 
two colorings to define two wavefunctions described by equation [101 |^co/)and l^^^)- 
We will use \Rcoi) and |-R ^) respectively as the reference configurations to define the 
transition loops. We have computed the loop distributions for both from Monte Carlo 
sampling of these wavefunctions. 

We can relate \^{z)) to the 0(1) model by considering \^coi)- In this case, since 
only the occupied links in Rcoi are colored red, the total loop length £ = 2 * Nb, and 
—Log{z) acts as an effective loop tension. Accordingly, we can relate 2; to in the 
0(1) model. As discussed in section [5l z = 1 sits in the dense phase, corresponding to 
the dimer liquid, and for some < 1 we expect to cross a phase transition into the 
dilute loop phase, corresponding to the dimer crystal. The relationship is less direct for 
l^v^)' ^tiere loop segments within the red hexagons are given unity weight, and only 
loop segments passing through black links are given a weight z < 1. However we still 
expect a transition from the long loop dense phase, to a short loop phase of finite loop 
density. 

In figure |5l the left panels show the total loop density p, and the right panels show 
the expectation value of the longest loop length, s^Im) as a function z, for |\l/co«)(top) 
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and ^) (bottom). Here we have rescaled s{iM) by L'^^'^ where 7/4 is the fractal 
dimension of the dimer hquid. In both cases s{iM) / L'^^'^ remains finite for z close to 1 
and vanishes for small z. Finite size scaling of s{£m) suggests that the corresponding 
phase transitions occur at z™' ~ 0.57 for l^^coi) and z/^ ~ 0.53 for |^'^). For I'ifcoi), 
p vanishes in the columnar phase, whereas p remains finite in the \/l2 x \/l2 phase 
of I^E'^/]^)- This distinction is due to the fact that in the columnar phase, dimers are 
pinned to \Rcoi), whereas the a/T2 x a/T2 phase dimers resonate within the hexagons and 
are not pinned to \R^)- The location of these transitions may be compared with the 
transition of the honeycomb lattice 0(1) model, where Kc = l/"\/3 ^ 0.58. Figure [6^ 
shows the longest loop distributions -Pfj,,^(s) for |\E' co«) (left) and |\E' ^) (right) for several 
values of z. We see that in both cases the broad distribution, which is a signature of 
the liquid phase, vanishes below the transition at Zc- The fractal dimension within the 
dimer liquid phase is plotted in figure Eb- Deviations from Df = 7/4 in the liquid phase 
near the transition may be due to the limitations of extrapolating the finite size scaling 
on lattices of these sizes. 

We note that in the crystalline phases, the loop distributions depend on the choice 
of l-Ro); here we have chosen reference dimerizations which reflect the known ordering 
in each phase. In particular, the vanishing of s{C.m)/L^^ in the crystalline phases is 
due to a choice of |-Ro) that is commensurate with the order. While in static crystals 
such as the columnar phase we may expect that any choice of |i?o) will lead to s{Im) 
scaling as L", L^, or L^, in a resonating dimer crystal such as a/T2 x a/T2, this is not in 
general true. In fact, we find that for an arbitrary choice of |i?o) the loop distribution in 
a resonating crystal phase may appear to be that of a liquid phase. This is illustrated 
in figure [7] where we have computed P^^^- for the ideal a/T2 x a/T2 state |^E'^(-2 = 0)) 
for two different columnar reference configurations: \Rcoi) (seen in figure [1^) and \Rcoi) 
which is related to \Rcoi) by translating all dimers one link horizontally. Figure [7] shows 
that for \Rcoi), Pi^ appears to be a liquid phase, but for|i?co«) it is sharply peaked at a 
single length scale. In fact finite size scaling of the length of Em shows that for \Rcoi), 
the fractal dimension is that of the liquid phase {Df = 7/4), but for \Rcoi) Df = 1; 
indicative of the crystal order. This suggests that one may distinguish a liquid phase by 
computing Df for all (in this case 12) configurations that are related by symmetry to 
l-Ro), to ensure that a resonating crystalline order is not hidden by the choice of |-Ro)- 
However, for any symmetry broken state, if \Ro) is chosen to be commensurate with the 
broken symmetry, the length of Im will be finite. 



7. Interacting dimer loop liquid 

Section [6] discusses the addition of an effective loop tension to the RK dimer liquid. 
Here we consider the effects of a translationally invariant dimer interaction by defining 
a generalization of the RK wavefunction: 

l*«) = ^E"~'^^^'^^l^) (12) 
VZ ^ 
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Figure 5. Top two panels: The total loop density p, and and longest loop length 
s(^A/) for l^coi)- Bottom two panels: The total loop density p, and and longest loop 
length s(^m) for l^^q^)- The solid lines are extrapolations to the thermodynamic 
limit. We have rescaled s{iM) by i^/'' for data collapse in the liquid phase. 




Figure 6. (a) Longest loop distribution P^j^- of j'l'coi) (left) and j^* (right) plotted 
for several values of z for a L = 60 lattice, (b) Fractal dimension -D/ of |^co;) and 
l^\/T2^ in the dimer liquid phase computed from the finite size scaling of s{£m) with 
lattice sizes up to L = 144. The line shows the value at the RK point, Df — 7/4. The 
accuracy of is limited by the lattice sizes studied. 
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0.2 0.4 0.6 0.8 1 

Figure 7. The longest loop distribution Pij^^ in the ideal \/T2 x \/T2 crystal, 
= 0)) computed for two different columnar reference configurations \Rcoi) 
and \Rcoi)- In a resonating crystal phase, for some choices of |i?o)j Pim may clearly 
indicate crystalline order (seen by the peak at a single length scale for the choice \Rcoi)) 
whereas for other choices (e.g. \Rcoi)), the distribution Pg^^, may appear to be that of 
a liquid. 



In f lT2|) . Nf (C) is the number of flippable plaquettes in |C). For a = 1, |$a=i) = \RK)', 
for a <^ 1, configurations with the maximum number of plaquettes are favored. The 
square norm of \^a) is the partition function of an interacting classical dimer model 

m- 

($J<I>«) =Z,, = ^e-'^-^-(^) (13) 

c 

Ea{C) =-uNf{C), « = exp(^^^. (14) 

All equal time correlation functions of \^a) are equal those of the classical model. This 
relationship allows us to write down a local RK-like Hamiltonian [l2l |26] for which |$q,) 
is an exact zero energy ground state: 

H^ = J2 + ^-c-) +t(a'^^f^^ |Z7)(Z7| + (15) 

p 

with ANf = Nf{/S/) - iV/(ZS7) and Nf{^) (iV/(A7)) is the total number of flippable 
plaquettes on the lattice with p in orientation /S7 (A^). Here we see that a = 1 ^ u = 0, 
which is the non-interacting point of the classical model; consequently Ha=i = -^rk at 
the RK point {v = t). The interacting classical dimer model has previously been studied 
using transfer matrix techniques by Trousselet et al. [13]. For u > (a < 1), Ed is 
minimized by configurations which have the maximum number of flippable plaquettes, 
Ni/6. There is a large number of such configurations, including the 12 symmetry related 
columnar states and all configurations related to these by translating rows of dimers 
(row shifting modes) or rotating all dimers in a set of columns (column shifting modes). 
For a = 0, this degeneracy prevents the formation of local order. However an "order- 
by-disorder" mechanism [Hj might allow for fluctuations to favor an ordered state for 
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nonzero a. Indeed the results of Ref. [43] are consistent with a first order phase 
transition from a hquid phase for a > etc to an ordered phase for a < ac, where the row 
shifting modes favor configurations with dimers ahgned in the same direction. 

Local updates do not give access to the global defects that appear as a — >■ 0. 
Therefore we have implemented a directed loop Monte Carlo algorithm ^51 HE] to sample 
|$Q,) on lattice sizes up to L = 128, using \Rcoi) as the reference configuration. The 
dimer model directed-loop algorithm generates non-local updates by creating a pair of 
defects that violate the hard core, fully packed dimer constraint; these defects undergo 
a directed random walk around the lattice until they coincide and annihilate [15]. The 
path of these defects forms a closed loop, along which the dimer configuration is updated. 
To sample |$q-), we have chosen the local weights of each step in the random walk to 



reflect |$q,) [47|, and found a solution to the directed-loop equations [46] (see Appendix A 



for details). This algorithm is ergodic in all topological sectors, and therefore the results 
presented here are computed over all winding sectors. 

Figure M shows that in the presence of interactions, the dimer liquid persists down 
to a ~ 0.2. For much smaller values of a, Dj approaches 1, indicative of a symmetry 
broken phase. However there is an intermediate regime where the finite size scaling of 
£m does not converge to a power law (this is shown by the dotted line in figure [8]). The 
longest loop distrbution P^j^^ (right panel of figure [8]) displays a series of peaks in this 
intermediate regime, whereas for small a, there is only a single peak at s(£m) = L. 
The small a regime is consistent with the low temperature phase of Ref. [43j where 
row shifting modes dominate; with this choice of reference dimerization row defects are 
loops of length L winding around the lattice in the direction of the rows. 

To further characterize the appearance of this symmetry breaking order for small 
a, we define a rotational symmetry breaking order parameter 

where i and j each represent one of the three possible dimer orientations on the triangular 
lattice, and is the number of dimers in a given orientation. Additionally, we define 
a columnar order parameter, 

<1 " ]^ E (^^ ^] - i^col'] ) ' , (17) 



d 



(r,c) 



where -R^^f^ is a columnar configuration with rows oriented in direction r and columns 
in direction c, and -R^of is obtained by translating all dimers by one lattice spacing in 
direction r. A^^^ [R] is the number of dimers that coincide with the configuration R. In 
a columnar ordered phase, the expectation values of both M^^^ and M^^j saturate to 1, 
whereas M^^^^ = 1 and M^^j = in the rotational symmetry broken ordered phase that 
is favored by the row shifting modes. 

Figure [9] shows the expectation values of these dimer order parameters in the 
symmetry broken phases. We see that for small a, M^^^ saturates to 1, whereas M^^y 
scales to zero for larger system sizes (figure [9l). This corresponds to the rotational 
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(a) (b) 

Figure 8. (a) Fractal dimension of the largest loop of \^a) computed with a 
directed loop algorithm for lattice sizes up to L = 128. The solid line corresponds to 
Df — 7/4, the value at the RK point. The dotted line is a guide for the eye; the gap 
in data is where the finite size scaling of s(£j\/) doesn't converge to a power law for 
the lattice sizes studied, (b) Longest loop distribution of on a L = 64 lattice, 
for several values of a, computed using \Rcoi)- For very small a, Pij^^ is peaked at L, 
consistent with an ordered phase with all dimers aligned in the same direction. 



symmetry broken phase where row defects dominate and destroy the columnar order, as 
discussed in Ref. [13]. However, there is an intermediate regime for larger system sizes 



where the expectation value of M'^^^^ reaches at plateau at approximately 1/4, and M^^j 
is finite and approaching 1/2. This is consistent with a phase for which column shifting 
modes dominate, such that two dimer directions are preferred. The ideal columnar 
defect configurations saturate the order parameter 



R 



(r,c) 



R 



(r',c) 



R 



(r,c) 



col 



R. 



(r',c) 



col 



:i8) 



J\j2 ^ ] y^^d y- ^col J ^ " d y ^col 
(r,r',c) 

to 1, while M^Q^. and M^^j are 1/4 and 1/2, respectively. In ( fTSll . the sum is over 
r ^ r' ^ c. Figure [3 shows that, in the intermediate phase, M^^ approaches 1 with 
decreasing a. We note that while perfect column defects cost zero (classical) energy, 
adding a horizontal kink to the column defect costs a finite energy; the increased 
degeneracy of kinked defects may favor them over perfect defects for sufficiently large 
a. If kinked column defects proliferate, they will generate a finite (but non-maximal) 
expectation value of M"^. This suggests that for a — )■ 0, global row defects dominate but 
there is an intermediate regime where kinked column defects dominate. Both of these 
transitions appear to be first order as shown by the discontinuities seen in these order 
parameters for larger system sizes (see figure [9]); we have confirmed this by observing 
a double peak structure in the histograms of these order parameters at the transitions. 
Finite size scaling of these transitions suggests that both transitions occur at finite a in 
the thermodynamic limit. 

The distinction between these results and what is seen in Ref. [13] is likely due to 
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the difference in the geometry of the systems studied. We have also studied triangular 
lattices with Ly, where L^^y are the lengths in the two lattice directions. We find 
that as Ly/Lx increases, the width of the intermediate regime decreases. This follows 
from the discussion of Ref. [13]: increasing Ly/L^ increases the number of row-shifting 
modes, and therefore favors the row defect ordering. Ref. [13] considers lattices with 
Ly ^ L^, so the results of that work are consistent with ours in this limit. 




0.1 0.2 0.1 0.2 0.1 0.2 0.3 



a a a 

Figure 9. Expectation values of dimer order parameters M^^^^ (left), M^Qj(center), 
M^jj(riglit) in \^a) in the symmetry broken phases. For larger values of a, all local 
order parameters vanish. 



8. d-isotopic quantum loop gas in the triangular lattice QDM 

So far we have only considered transitions out of the 0(1) dimer loop liquid to crystalline 
phases; we now consider a wavefunction that allows us to explore other loop gases in the 
triangular lattice QDM. In reference [12] Freedman et. al. introduced a generalization 
of the quantum loop gas of the toric code they termed the (i-isotopic loop gas [l2]. By 
giving a weight d to each contractible loop in the ground state wavefunction, the 0(1) 
loop liquid of the toric code at d = 1 can be driven to loop gas phases for 0{n > 1): 

c 




Figure 10. The decorated triangular lattice from reference [12] used to define 
Here the staggered configuration defined by the red links is used as the reference 
dimerization. 
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In equation f ll9p . A/'c is the number of contractible loops in \C) and the sum over C is 
taken over all configurations in the same winding sector. The norm of is equal to 
classical partition function of the 0{n) loop model for n = (P [IE]. Freedman et. al. 
also wrote down an RK-like Hamiltonian, for which 1^^) is the exact zero energy 
ground state [12]. Every winding sector contains a unique ground state; as such this 
d-isotopic loop gas has a large degeneracy. For 1 < < 72, l^rf) is a critical loop 
liquid and it has been shown that is gapless. For n > 2, the O (n) model is in a 
short looped phase [SS], and correspondingly, the d-isotopic loop gas is in a short looped 
phase for d > \/2. All diagonal local correlation functions in vanish exponentially; 
however Troyer et. al. showed certain that off-diagonal correlation functions have a 
power-law decay [19]. The original motivation for introducing 1"^^) was the search for 
gapped topological phases beyond the abehan toric code phase. Subsequent work [49] 
shows that adding local interactions to Ha to open a gap drives the loop gas to the toric 
code phase. Consistent with the disordered nature of the liquid phase of Troyer 
et. al. showed that the critical phase of the loop liquid could be described by the fractal 
dimension of the long fractal loop |49j . 

To implement in a quantum dimer model we color the triangular lattice 
according to reference [12]. Links occupied in a staggered reference configuration \Rq) 
are colored red and links connecting the red links are colored green (see figure [10] ). 
Here we choose a slightly different interpretation of the transition loop graph; dimers 
in |C) that coincide with \Rq) are considered to be minimal, length 2 loops. As such, 
the loop coverings defined by the transition graphs are fully packed loop coverings. In 
this case, is described by a fully packed loop model, which in turn is related to the 
dense phase of the 0{n) model [50]. Consequently, the geometric exponents as defined 
in section [5] are the same in the fully packed loop model as those of the dense phase of 
the 0{n) model with d^ = n [10| |5T]. 

With this coloring, a local quantum dimer model Hamiltonian, Hd can be defined 
such that \^d) defines the unique, zero energy ground state in each winding sector. 
To allow for the transition to the short looped, staggered crystal phase in the limit 
(i — )■ oo, 4-dimer dynamics that fiip dimers around a length 8 loop are included in Hd 
to connect the staggered configurations to other configurations. Additionally, following 
reference [12], we wish to define Hd such that it only connects configurations in the same 
loop winding sector, not in a parity defined winding sector. This requires removing the 
plaquette fiip terms that can act as a surgery and connect different winding sectors-this 
is done by excluding plaquette flips in Hd for rhombi that have a green diagonal link. 
To compensate for a lost ergodicity due to the removal of these plaquette flip terms, a 
3-dimer resonance around a length 6 triangular loops is also added to Hd [12] . 

Figure [Tib shows the longest loop distribution for several values of d computed by 
Monte Carlo sampling of \^d) in the triangular lattice quantum dimer model. We see 
that a liquid phase persists for d > 1. Figure [TTk shows that the fractal dimension is a 
continuously varying function of d. Thus 1$^) can describe dimer liquid phases distinct 
from the Z2 topologically ordered phase at the RK point. 
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(a) (b) 

Figure 11. (a) Longest loop distribution Pij^,, for \'^d) in the zero winding sector on 
a L = 24 lattice for several values of d. (b) Fractal loop dimension of \'^d), computed 
from finite size scaling of s{£m) on lattice sizes up to L = 84. 

9. Off-diagonal loop operators in the triangular lattice QDM 

Off-diagonal operators in a quantum dimer model that do not violate the hard-core 
constraint are loop operators: two dimer configurations can be connected by flipping 
dimers along their transition loops, as given in equation ([6]). In the toric code, such 
loop operators are related to the Wilson loop of the underlying Z2 gauge theory. In, 
the ground state |^tc), the expectation value of all loop operators is exactly 1 (a "zero 
law"). However the addition of a magnetic field will cause the loop operators to decay 
with the length of the loop (a "perimeter law"). In a quantum dimer model, the fully- 
packed constraint means that {Wi{s)) < 1 for all loops, even at the RK point [2^. 
Moessner, Sondhi and Fradkin suggested that (PV^(s)) will decay as a perimeter law at 
the RK point due to the extensive entropy of dimerizations of the lattice [28] . Figure [T2l 
shows (VF^)rk computed for a set of parallelograms of differing lengths on the triangular 
lattice by Monte Carlo sampling of \RK). {We{s)) is seen to display a clear perimeter 
law at the RK point. 

10. Conclusions 

In this paper we have characterized the formation of topological order in quantum dimer 
models in terms of the physical mechanism of loop condensation. This allows for a more 
direct connection between the Z2 topologically ordered liquid of the triangular lattice 
QDM to that of the toric code and other lattice models. We have shown that geometric 
properties of the transition graph loop distribution may be used to distinguish dimer 
liquid and crystal phases. Additionally, the fractal dimension of the spanning loop may 
be used to distinguish dimer liquid phases. We believe this approach may be useful for 
the study of dimer liquids and the full phase diagram of the canonical quantum dimer 
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Figure 12. Expectation value of off-diagonal loop operator We as a function of loop 
Icnth s at the RK point, computed for various parallelograms of different lengths on 
lattices of size L — 16, 32. The line shows a fit to an exponential decay, indicating a 
clear perimeter law behavior. 

model Hamiltonian. 
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Appendix A. The directed-loop algorithm in the interacting dimer model 

The directed-loop algorithm [TTl H5] generates non-local updates in a dimer model by 
creating a pair of defects (a doubly occupied vertex and an unoccupied vertex) that 
violate the fully-packed, hard-core dimer constraint. One of these defects undergoes 
a directed random walk and, when the two defects coincide once again, this generates 
a new allowable dimerization, where the dimers have been flipped along the path of 
the path of the defect |l5l SB]. To maintain detailed balance, the probabilities of 
choosing each direction on a given step of the random walk must satisfy the directed-loop 
equations |i7] . 

Consider a step in the loop-update where the loop has entered a vertex v from link 
/. We define the probability of choosing link /' to exit the vertex as Pi^u = ti^i'/wi, where 
Wi is the weight of the configuration with the dimer on link /. Detailed balance and 
probability conservation require 

^ti^ii = wi, ti^v = tv^i (A.l) 
/' 

For a given set of weights {w^}, this leads to a set of directed-loop equations. 

Here we consider sampling the square of the wavefuction = a~^fp. There are 

three possible weights for a dimer to occupy a link {wq, Wi, 11)2}^ corresponding to a link 
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with 0, 1 or 2 parallel dimers. The choice W2 = 1, so W2 = wf = gives the following 
directed-loop equations: 

1 = (6- ni-n2)t02 + r^ltl2 + (^2-1)^22 + ^2 (^-2) 

^2 = (6 - m - na) toi + n2ti2 + K - 1) tu + (A.3) 
= ^^^^^ _^ ^2to2 + (6 - ni - ria) too + ^q- (A.4) 

In (1X211X41) . is the weight for a transition between links with p and q parallel dimers, 
Up is the number of links with p parallel dimers at vertex v, and is rate to "bounce" 
back from a link with p parallel dimers. Since these equations are underdetermined, we 
can choose a solution to flA.2tiA.4| for which the bounce probabilities vanish at the RK 
point (a = 1): 



too — toi — to2 — ti2 — — , tg — 0, (A. 5) 



til = { — ^1 > 1 , (A.6) 

ni < 1 



a'^n2-a^il-a^ ni+6 1-a"' -1 



''22 



^2 > 1 



sR^y "2 - " , (A.7) 

n2<l 

( «^(l-«^)(5-na) _ 

1 n,^l 

tl f l-*(5«^ + (l-«>i) n2 = l _ (^g^ 
I n2^l 

The directed loop algorithm using this solution to (lA.2tlA.4| has be compared to a local 
update algorithm and is sufficiently efficient in the regimes studied. 
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